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1 Introduction 



This paper continues investigation of our previous paper 1 where we studied 
the dynamics of the cosmological model, including gravity and minimally coupled 
scalar field with simple potential including only massive term. The dynamics of the 
minisuperspace cosmological models with the massive real scalar field for the flat, 
open and closed Friedmann universes was studied in papers 2 ' 3 in terms of phase 
space. It appeared that the dynamics of closed model is more complicated than 
that of open and flat models. This dynamics allows the transitions from expan- 
sion to contraction and the existence of points of maximal expansion and minimal 
contraction in contrast with cases of open and flat cosmologies. Moreover, closed 
spherically symmetric models cannot expand infinitely and should have the points of 
maximal expansion provided the matter in the model under consideration satisfies 
the condition of energodominance 4 . The presence of points of maximal expansion 
and minimal contraction opens the possibility for the existence of the trajectories 
escaping singularity and oscillating between turning points 5 . The possibility of 
existence of such trajectories or non-singular universes filled with scalar field was 
discussed also earlier in Ref. 6. Recently such a dynamics was studed from the 
viewpoint of theory of dynamical chaos 12 . 

Our previous work 1 was devoted to the classification of trajectories in closed 
model and to the study the idea about the fractal set of infinitely bouncing uni- 
verses which was put forward in Ref. 5. Here we consider more rich model where 
cosmological constant is included into consideration. Dynamics in such a model is 
rather peculiar and is of interest from mathematical point of view. Besides, some 
attempts to reconstruct the potential of inflaton field studying the anysotropy of 
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relic cosmic background radiation testify in favor of the existence of constant term 
in the power expansion of this potential 7 : 

Apparently, V plays role of cosmological constant A. 
Thus, we have an action 

S = J d^x^—g^{R - 2A) + l -g^d^d vV - in* 2 ? 2 } , (1-1) 

where rap is Planck mass. The presence of cosmological constant A makes possible 
an existence of trajectories without points of maximal expansion, i.e. the trajecto- 
ries which begin and end in DeSitter regime 

a(t) <~ exp(— Ht), t — > — oo; 

a(t) ~ exp(ifi), t -> oo (1.2) 

or the trajectories which begin in singularity and end in DeSitter expansion or vice 
versa. Here, in Eq. (1.2) H denotes Hubble constant 




It is convenient to begin the consideration of possible classes of trajectories with 
an investigation of the simplest cosmological model when scalar field is massless, 
but cosmological constant A ^ (the model with m = 0, A = was studied in 
some detail in Refs. 8,1,9). The second section of our paper will be devoted to this 
case. The class of trajectories travelling between sigularities with possible multiple 
bounces will be described in third section in terms of configuration space. In fourth 
section we consider the trajectories of two types DeSitter-singularity and DeSitter- 
DeSitter and the surface separating them in phase space. 
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2 Massless model with cosmological constant 



Equations of motion can be now written down in the following form 

y = -3yz, (2.1) 

and 

, A 8-Ky 2 

Z = Z+ 3-^ (2 ' 2) 

where we have introduced notations 



a 

y = <p, z = - 



and dots denote differentiation in respect with time. First integral of motion looks 
as follows: 



1 _ 47TJ/ 2 A 



a 1 3m z p 3 



As far as the variable <p is not included into Eqs. (2.3)-(2.5) the corresponding 
dynamical system can be represented by the diagram of two-dimensional phase 
space with coordinates (y, z), see Fig. 1. Here two hyperbola 



separate the part of phase space accessible for the trajectories of closed and open 
Friedmann - Robertson - Walker model (cf. Refs. 2,3). There are four special 
points: points A and C are attracting and repelling nodes respectively, while points 
B and D are saddle ones. The coordinates of these points have the following values: 



y A = 0, z A = y -; ys = y -g^j - ! z b = 0; 



/A m 2 p A 
yc = 0, z c = -]j-; y D = -y-^-, z D = 0. (2.5) 



The region I bounded by trajectories connecting the points A and B, B and C, 
C and D, D and A is covered by trajectories which begins at point C, i.e since 
contraction according to DeSitter exponential law at t — — oo, then have a bounce 
or in other words go through the point of minimal contraction and end in the point 
A corresponding to DeSitter expansion. 

The region 77 contains trajectories expanding from singularity until DeSitter 
regime while the region III contains trajectories contracting from DeSitter point C 
to singularity. Finally, the region IV contains trajectories which begin expansion 
from singularity, then go through the points of maximal expansion and then return 
to singularity again. The regions II and III are separated from the region IV by 
trajectories connecting saddle points B and D with singularity. 

It is interesting to notice that the picture analogous to that described above 
was obtained earlier for the hydrodynamical matter in the presence of cosmological 
constant in Ref. 10. The diagram similar to diagram in Fig.l was constructed in 
the coordinates e (energy) and z. The stationary unstable solutions were obtained 
in Ref. 10 as well. Such a coincidence is quite natural because massless scalar field 
corresponds to the hydrodynamical matter with equation of state e = p, where p is 
pressure. 

Inclusion of mass of scalar field makes the picture more complicated. Phase 
space becomes 3-dimcnsional and the opportunity of multiple bounces or passages 
via points of maximal expansion arises. Besides, the saddle points B and D disap- 
pear. Nevertheless, one can speak roughly about three regimes again: they can be 
called "singularity-singularity" , "singularity-DeSitter" and "DeSitter-DeSitter" . 
First type of trajectories correspond to the trajectories populating region IV in 
massless case, but their behaviour can be more involved due to passage through 
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points of minimal contraction and maximal expansion. Classification of these tra- 
jectories in terms of configuration space (a, ip) will be given in third section of this 
paper. The second type of trajectories corresponds to those in regions II and III 
in the massless case while "DeSitter-DeSitter" trajectories correspond to those from 
the region I. These two types of trajectories will be studied in section 4. 

3 Travelling between singularities 

It was shown in the preceding paper 1 that the points of maximal expansion 
and those of minimal contraction can exist only in Euclidean region. In the case of 
the presence of cosmological constant A the form of this region is given by equation 

2 3m 2 P Amp 
V < 1^~ 4^" (3 - 1} 

and this region is closed on the right at 

The area where located the possible points of maximal expansion is separated from 
that with points of minimal contraction by the curve 

, _ m% _ Arn^ 

P — n 2 2 A 2 ' V J - Z J 

which crosses the horizontal axis at 

= VT 

Thus, geometry of this curve surrounding the region of possible points of maximal 
contraction differs from that consisting of two branches in the absence of A (cf. Ref. 
f). The configuration of Euclidean region and the separating curve ( |3.2|) is shown 
in Fig. 2. 
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In the preceding paper 1 it was shown that the region of possible points of 
maximal expansion has quite a regular structure. In the left, closest to axis a = 
part of this region there are points of maximal expansion after which trajectory goes 
to singularity. Then one can see the region where after the going through the point 
of maximal expansion a trajectory undergoes the "bounce" i.e. goes through the 
point of minimal contraction. Then we have the region where after going through 
the point of maximal expansion trajectory has a "<p - turn " i.e. has the extremum 
in the value of the scalar field ip and then fells into singularity. Then one has the 
region corresponding to the trajectories having bounce after one oscillation in <p and 
so on. Then studying the "substructure" of regions corresponding to trajectories 
possessing bounces, one can see that this substructure repeats the structure itself. 
In turn the "subsubstrucure" of regions having two bounces repeats the general 
structure and so on and so forth. Continuing this process ad infinitum we can get 
the set of infinitely bouncing trajectories escaping the singularity. It is plausible to 
believe that this set has a fractal nature 11 due to recurrent type of the procedure 
of its construction s* 1 ' 9 * 12 . 

The remarkable feature of the model with the cosmological constant A consists 
in the fact that this fractal structure of the set of infinitely bouncing trajectories 
survives provided the value of A is small enough in comparison with the mass of 
scalar field. 

An important distinction between the case of small A ^ and A = consists 
in the fact that an area of all the bouncing intervals (i.e. such locations of points of 
maximal expansion after which the trajectories have bounce) is restricted from the 
right in configuration space by the new unstable periodical trajectory which has no 
analogue in the case of A = 0. Trajectories situated on the left side of this periodical 
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trajectory go to the regions of possible bounces while those situated on the right go 
to DcSitter point. Another distinction is that substructure of all intervals contains 
subintervals, corresponding to trajectories, going to DcSitter after bounce. 
At some critical value of 

A - 0.5m 2 

given by the numerical computation the complicated structure of subregions of 
the part of Euclidean region containing the points of maximal expansion vanishes. 
Instead one has only two regions (see Fig. 3): the left one containing the points 
of maximal expansion after which trajectories fell to singularity and the right one 
where after going through point of maximal expansion and then via point of minimal 
contraction trajectories go to DeSitter regime. Thus for the case of large values of 
A the behaviour of trajectories qualitatively reminds that for the case of massless 
scalar field described in the preceding section. 

It is interesting to notice that the transition from the chaotic (fractal) structure 
of trajectories when we have infinite number of subregions to the regular structure 
(Fig. 3) with two subregions has jump-like character. There is not transitional 
situation when one can see finite number of different subregions. 

4 Configuration of trajectories in the phase space 

In this section we consider the relation between the trajectories travelling be- 
tween two DeSitter regimes and those which begin at DeSitter point and end in 
the singularity or viceversa. The boundary between these two classes of trajectories 
is rather complicated and fractal-like two dimensional surface in three-dimensional 
phase space. To describe it we shall use smooth two-dimensional projections ignor- 
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ing its fractal nature. 

It is convenient to introduce the following coordinates in the phase space: 

mpip a 

x = mpip, y = , z = m— . 

to ' a 

At A 7^ we have two De-Sitter points at x = 
V = « 

z = im-^/A/3 Some part of trajectories beginning in the neighborhood of DcSittcr 
points go to another DeSittcr point while the rest of trajectories fell to singularity 
or in terms of phase space to an infinite sphere (see Refs. 2,3). To distinguish 
between these two classes of trajectories we make numerical integration of equations 
of motion starting from lower hemisphere around De-Sitter point. For concreteness, 
the radius of this hemisphere is chosen as r — z/3. The behaviour of trajectories is 
shown in Fig. 4. It is convenient to parametrize initial conditions for trajectories 
on hemisphere by two angles <p an d which are defined as 

tan(</>) = y/x; sin(0) = (y/A/3 + z)/r. 

While f^ooa trajectory may go to singularity or to DeSitter asymptotic point. In 
Fig. 4 dashed area corresponds to trajectories going to singularity. It is necessary to 
emphasizes that this picture is a rough one - it exists a fine structure of trajectories 
escaping sungularity inside the dashed area due to bounces. In Fig. 4b such a 
structure is remarkable. For every value of relation A/to 2 there are two distinguished 
values of an angle 9: Q m in such that at 9 < 9 m in all the trajectories go to singularity 
independently of value of angle <f> and 9 max such that all the trajectories go to the 
second DeSitter point while 9 > 9 max . 

It is interesting to trace out the dependence of these angles on the relation 
A/to 2 . At small values of A/m 2 these angles almost coincide 9 m in ~ 9 max- This 
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(4.1) 



fact can be explained by the fact that the distance between DeSitter points in 
the phase space is large in comparison with mass parameter and the trajectories 
manage to have so many rotations around axis z that the dependence of their final 
destiny on the initial angle 4> disappear. It would be reasonable to expect that in 
the limit A/m 2 — ► oo these two limiting angles should also coincide giving the angle 
separating different type of trajectories in massless case (cf. sec. 2), however it is 
not the case because in our parametrization small values of mass m are compensated 
by possible large values of the field tp and one does not have smooth transition to 
the massless model. Therefore, it is interesting to carry out numerical investigation 
of massless model considering since the beginning two-dimensional phase space. In 
this case one has only one parameter for the trajectories going from lower DeSitter 
point. It is "latitude" angle 9. The value of 9 separating trajectories going to 
singularity in the limit r — > appears to be equal 



' n umcr /cat 



(4.2) 



One can easily estimate theoretical value of this angle connecting DeSitter point 
with one of two saddle points of two-dimensional phase diagram (see Fig. 1). Using 



Eqs. (2.5) together with normalization given by Eqs. (4.1) one can easily get 



Otheor = arctanW w 1.3. (4.3) 

Thus we have seen that the results of theoretical consideration and numerical sim- 
ulation are in a good agreement. 
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Captions to Figures 
Fig. 1 In Fig. 1 the phase diagram for the model with massless scalar field is 
presented. Two hyperbolah going thgrough points A and C separate the part of 
phase space accessible for the trajectories of closed and open Friedmann - Robertson 
- Walker model. Points A and C are attracting and repelling nodes respectively, 
while points B and D are saddle ones. The region / bounded by trajectories con- 
necting the points A and B, B and C, C and D, D and A is covered by trajectories 
which begins at point C, i.e since contraction according to DeSitter exponential 
law, then have a bounce and end in the point A corresponding to DeSitter expan- 
sion. The region II contains trajectories expanding from singularity until DeSitter 
regime while the region III contains trajectories contracting from DeSitter point C 
to singularity. Finally, the region IV contains trajectories which begin expansion 
from singularity, then go through the points of maximal expansion and then return 
to singularity again. The regions II and III are separated from the region IV by 
trajectories connecting saddle points B and D with singularity. 
Fig. 2. Solid line denotes the boundary between Lorentz and Euclidean regions, 
dashed line corresponds to the curve separating the points of minimal contraction 
and maximal expansion. 

Fig. 3. Separating curve (solid) and boundary of bouncing interval (dashed) in 
the case of large A/to 2 . The bouncing interval is located between these two curves. 
Fig. 4 In Fig. 4 dashed area corresponds to the trajectories going out DeSitter 
contracting point and falling to singularity. The angles 9 and ip parametrize the 
direction of the trajectory in the initial DeSitter point. White area corresponds to 
trajectories going between two DeSitter points. In Fig. 4a A = 0.05m 2 , in Fig. 4b 
A = 0.5m 2 , in Fig. 4c A = 50m 2 . 
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